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Abstract 

In this paper we study the singularity of multivariate Hermite interpolation 
of type total degree. We present a method to judge the singularity of the inter- 
polation scheme considered and by the method to be developed, we show that all 
Hermite interpolation of type total degree on m = d + k points in W d is singular 
if d > 2k. And then we solve the Hermite interpolation problem on m < d + 3 
. nodes completely. Precisely, all Hermite interpolations of type total degree on 

\ m < d + 1 points with d > 2 are singular; for m = d + 2 and m = d + 3, only 

three cases and one case can produce regular Hermite interpolation schemes, re- 
spectively. Besides, we also present a method to compute the interpolation space 
\ for Hermite interpolation of type total degree. 



Keywords: Hermite interpolation; Singularity; Interpolation space; Polynomial 
ideal 



>< 

$ : 1 Introduction 

Let Tl d be the space of all polynomials in d variables, and let il^ be the subspace of 
polynomials of total degree at most n. Let = X 2 , ■ ■ ■ , X m } be a set of pairwise 
distinct points in M d and p = {p±,P2, ■ ■ ■ ,p m } be a set of m nonnegative integers. The 
Hermite interpolation problem to be considered in this paper is described as follows: 
Find a (unique) polynomial / 6 11^ satisfying 

Qa 1 +a2+...+a d 

dxT ...dx ad f< " Xq) = C<? ' Q ' 1 - q ~ m ' - |c *' - Pq ' (1) 



*This project is supported by NNSFC (Nos. 61033012,11171052,11271060,61272371) and "the Fun- 
damental Research Funds for the Central Universities" . 
^Corresponding author: mzhl@ dlut.edu.cn 
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for given values c ?)Cn where the numbers p q and n are assumed to satisfy 



Following [11|, |12j, such kind of problem is called Hermite interpolation of type 
total degree. The interpolation problem (p, is called regular if the above equation 
has a unique solution for each choice of values {c qiCl , 1 < q < m, < \et\ < p q }. 
Otherwise, the interpolation problem is singular. As shown in Q, the regularity of 
Hermite interpolation problem (p, JT) implies that it is regular for almost S£ C M d 
with \ 3£\ =m. 

Definition 1 (0). We say that the interpolation scheme p is: 

• Regular if the problem (p, is regular for all SC . 

• Almost regular if the problem (p, 3£) is regular for almost all 3£ . 

• Singular if (p, S£) is singular for all !% ' . 

The special case in which the p q are all the same is called uniform Hermite interpo- 
lation of type total degree. In the case of uniform Hermite interpolation of type total 
degree, Eq. (j2J) should be changed to 

n + d\ fp + d\ 

ml . (3) 



d J \ d 

The research of regularity of multivariate Hermite interpolation is more difficult that 
Lagrange case, although the latter is also difficult. One of the main reasons is that Eq. 
(J2J) or (J3]) do not hold in some cases. Up to now, we have known that all the Hermite 
interpolation on m < d + 1 points are singular except for Lagrange interpolation, see 
@, EH Qjj]- Besides, no any other results appeared for m > d + 2. Actually, Hermite 
interpolation of type total degree on d + 2 nodes in M. d are not necessary singular. For 
more research of this area, we can refer to [l|, 0-0, [8~13] and the reference therein. 



The main purpose of this paper is to investigate the singularity of Hermite interpo- 
lation for m = d + k with k = 1, 2, 3. Our method is to construct a polynomial which 
is a solution of the homogenous interpolation problem. By the presented method, we 
show that all Hermite interpolation of type total degree on m < d + 1 nodes in R d 
are singular except for Lagrange interpolation; on m = d + 2 nodes in W 1 are singular 
except for three cases; on m = d + 3 nodes are singular except for one case. The result 
of m < d+ 1 is well known, but our method is different. Moreover, we also show all the 
hermite interpolation problem of type total degree with m = d + k nodes are singular 
for d > 2k. 

This paper is organized as follows. In section 2, we will consider the interpolation 
space satisfying the Hermite interpolation requirement from the view of polynomial 
ideal. In this section, Eq. (j2J) is not required and the polynomial space is not necessary 
n^. In section 3, we consider the singularity of the Hermite interpolation of type total 
degree and present the main results. Finally, in section 4, we conclude our results. 



2 



2 Interpolation Space 



In this section, we consider the following interpolation problem: 
Let S£ = {Xi, X 2 , . . . , X m } be a set of pairwise distinct points in M d and p = 
{V11V21 ■ ■ ■ ,Pm} be a set of m nonnegative integers. Find a subspace V C IT d such 
that for any given real numbers c q>a , 1 < q < m, 1 < \ot\ < p q there exists a unique 
polynomial f EV satisfying the interpolation conditions 

Qai+ot2+-..+a d 

dxT...dx a / f{Xq) = Cq ^ l - q ~ m > °^\ a \^Pi ( 4 ) 

where ex = (ai_ L a 2 , • • • , &d) G Nq and \et\ = a\ + a 2 + . . . + Qd- 

Following [14| . we call such a pair {%~,p,V} correct. Clearly, such kind of space 
V always exists if no any constraint is added. For the Lagrange case, such kind of 
interpolation problem was studied extensively by many authors, for example ji], 0, @, @] 



and the reference therein. In [14|, Xu presented a solution from the view of polynomial 
ideal. Now we generalize Xu's result to Hermite case. 

Before proceeding, we first present some necessary notations. Throughout of this 
paper, we use the usual multi-index notation. To order the monomials in Il d = 
WL[x\, . . . ,xa], we use graded lexicographic order. Let J be a polynomial ideal in Il d . 
The codimension of I is denoted by codimi", that is, 

codimJ = dim U d / 1. 

If there are polynomials /i, f 2 , ■ ■ ■ , f r such that every / 6 I can be written as 

/ = ai/i + a 2 f 2 + • • • + a r f r , aj G Il d , 

we say that / is generated by the basis fi, f 2 , . . . , f r , and we write / = (fi, f 2 , . . . , f r ). 

For a fixed monomial order, we denote by LT(f) the leading monomial term for any 
polynomial / G H d ; that is, if / = c a X a , then LT(f) = &X^ , where X^ is the leading 
monomial among all monomials appearing in /. For an ideal / in U d other than {0}, 
we denote by LT(I) the leading terms of /, that is, 

LT(I) = {cX a |there exists f el with LT(f) = cX a }. 

We further denote by (LT(I)) the ideal generated by the leading terms of LT(f) for 

all/e/\{()}. 

The following theorem is important for out purpose. 



Theorem 1 (|l4|). Fix a monomial ordering on U d and let I e U d be an ideal. Then 
there is an isometry between U d /I and the space 

5/ := Span{X a |X a £ (LT{I))}. 

More precisely, every f G Tl d is congruent modulo I to a unique polynomial r G <S/. 
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Let / = (fx, / 2 , . . . , f r ) and J = g 2 , . . . , g s ) be two polynomial ideals. The sum 
of I and J, denoted by I + J, is the set of / + g where / G I and g G J. The product of 
J and J, denoted by J • J, is defined to be the ideal generated by all polynomials / • g 
where / G I and g G J. It is easy to know that I ■ J = (figj : 1 < z < r, 1 < j < s). The 
intersection J fl J of two ideals / and J in U d is the set of polynomials which belong 
to both I and J. We always have / • J C / fl J. However, IJ can be strictly contained 
in / fl J. It follows from [2| that if I and J is comaximal, then I J = I D J. / and J is 
comaximal if and only if / + J = n d . 

In application, people usually are interested in the space with minimal degree for 
fixed monomial ordering. For this purpose, consider the following polynomial ideal: 

/(I,p) = {/etf:- __/pg = 0, l<q<m, < \ct\ < p q } (5) 

If only one point Ifl" then we can write it as I(X,p). 

Theorem 2. Let I(JT,p) be the polynomial ideal defined as above. Then the interpo- 
lation problem satisfying Eq. (J4]) has a unique solution in 

Proof. Let N = YZ=i i^)- Denote the N linear functional in fl3J) by F l7 F 2 , F N . 
Thus Eq. (TjJ can be rewritten as F q (f) = c q for 1 < q < N. Suppose n is an integer 
which is big enough such that there exists a polynomial / G 11^ satisfying F q (f) = c q . 
We also assume that Si C LT^ although it is not necessary for our proof. Thus we have 
11^ = Si U (LT^ fl /). Suppose (fx, . . . ,ip t and 0i, . . . , S are the basis functions of Si and 
R I, respectively. We further define a column vector 

$ = (<fx, . . . , (p t , 0i, . . . , <p s ) T := ($f , $£) T . 

Since n is big enough, we have 

rank(F 1 ($),F 2 ($),...,F JV ($)) = JV. 

Furtherly, 

rank^^), F 2 ($i), . . . , Fjv($i)) = iV 

since Fj($ 2 ) = for 1 < i < N, which leads to t > N. It only remains to prove t < N. 
If t > N, then T± := (F 1 (<p i ), F 2 (<fi), . . . , F N (cpi)), 1 <i <t are linearly dependent and 
there exist scalars a±, a 2 , . . . , a t , not all zero, such that Yll=i a i-^i = 0- I n terms of the 
components of the vector this shows that 

t 

'£a i F q (<p i ) = Q, q = l,2,...,N 

i=i 

or, 

i 

F g (^a^)=0, q = l,2,...,N. 

i=i 

The latter equations means that <p = Ylt=i a i L Pi ^ Fy^,p), which is a contradiction 
to ip G Si because every ^ G 5/. Hence we have t < N and finally t = N, which 
completes the proof. □ 
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The theorem states that (i£*,p, <S/) is correct. This result maybe was known more 
or less, but we did not find it in the literature. 

Next, we consider the computation of I(3£ , p). If only one point is in 3£ , the result 
is immediate. 

Lemma 1. Let X be a point in M. d and p be a nonnegative integer, then 

I(X,p) = (l\h ■ ■ ■ lp+i : k is a linear polynomial vanishing at point X^. (6) 

For multi-point case, we also have the similar result. 
Theorem 3. 

/( JT, p) = ( / E n d :/ can be divided by the product of pi + 1 

(7) 

linear polynomials which pass through X i} i — 1, 2, . . . , m 



Proof. Without loss of generality, we only give a proof for m — 2, that is, X = {X±, X2} 
and p = {pi,p 2 }- Let 

I(X 1 , Pl ) = [f e U d : — "q^tJ{Xi) = 0, < |a| < Px }, 

I(X 2 ,p 2 ) = {/etf: — dxT f(X 2 ) = 0, < |a| < p 2 \, 

Obviously /( JT, p) = I (X ly pi) D I (X 2 , p 2 ) . If we denote by J the right hand of Eq. ([7j), 
then we need to show I(X 1 ,pi) fl I(X 2 ,p 2 ) = J. Obviously J C I(Xi,pi) D I(X 2 ,p 2 ). 
Thus it remains to show J D I(Xx,pi) fl I(X 2 ,p 2 ). Notice that I(Xi,pi)I(X 2 ,p 2 ) C J 
holds. Hence the proof will be completed if we can show 

i{x x , Pl ) n i(x 2 , P2 ) = i{x 1)Pl )i{x 2) p 2 ). 

To this end, we will prove that I(Xi,pi) and I(X 2 ,p 2 ) are comaximal, that is, 
I(X uPl ) + I(X 2 ,p 2 ) = IT*. It is enough to show 1 E I(X uPl ) + I(X 2 ,p 2 ). 

Let Xi = (x];, xjj,, . • • , x^) and X 2 = (x\ , x|, . . . , x^). Assume x{ 7^ xf. It is easy to 
check that (xi - x^) pi+1 E J(Xi,pi) and (x x - x?) P2+1 E I(X,p 2 ). If (x x - x\) Pl+1 and 
(xx — xf) P2+1 are seen as two polynomials with respect to Xi, that is, they are taken as 
univariate polynomials, then the greatest common divisor 

gcd((x 1 -xir +1 ,(x 1 -x?r+ 1 ) = l, 

which means that there exist two polynomial qi{xi) and ^2(^1) such that 

(xx - x\)^+ 1 q 1 {x 1 ) + {xi - xlr +l q 2 { Xl ) = 1. 

This complete the proof. □ 
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The following two examples will be mentioned again in next section to show the 
regularity of Hermite interpolation problem. 

Example 1. Consider the case of d — 3, m — 5 and p = {1, 1, 1, 1, 1}. Take 

X 1 = (0,0,0), X 2 = (1,0,0), X 3 = (0,1,0), X 4 = (0,0,1), X 5 = (x ,y ,zo) 

and «£" = {Xi, X 2 , X 3 , X 4 , X 5 }. With the help of Maple, the Groebner basis of /(#", p) 
can be written as 

{/(x, y, z, x , y , zo), f(z, y, x, z , y , x ),f(x, z, y, x , z , y ), 
g(x, y, z, x , yo, z ),g(z, y, x, Zq, y , x ),g(x, z, y, x , z , y ), 
g(y, x, z, ?/o, x , z ),g(z, x, y, z , x , y ),g(y, z, x, y , z , x ), 
h(x, y, z, x , j/o, z ), h(y, x, z, y , x , z ), h(y, z, x, y , z , x ), 

w(x, y, z, x , 2/o, zo),w(y, z, x, 2/0, zq, x q ),w(x, z, y, x , z , 2/0)}, 

where 

f(x,y,z,x ,y ,z ) = -x z (z - l)(x z + l(X 5 ) -y )(zy 2 + z 2 y - zy) 
+2(z - l)z (y z 2 + x z 2 - y z + x y z - x z + x y )xyz 
-y (z - 1)^0(2/0^0 + ^5) - x )(xz 2 + x 2 z - xz) 
+x y l(X 5 )(z 4 - 2z 3 + z 2 ) - z 2 (z - l) 3 (x 2 y + xy 2 - xy), 

g(x,y,z,x ,y ,z ) = x zo(x zo + z - l)(zy - zy 2 ) + z 2 (z - l) 2 (xy-xy 2 - x 2 y) 

+z (-z 2 + 2y z 2 + 2x0^ + 2x oVoZo + 2z - 3y z - Ax z - 1 + y + 2x )xyz 
+x l(X 5 )z 3 y - x (x z 2 + x Q + 2/0 + z 2 - l)z 2 y + (xz - x 2 z - xz 2 )y 2 z 2 

h(x, y, z, xq, 2/0, z ) = 2y z (y z - y + x y - z + x z + 1 - 2x )xyz 
-yoZo(z - l)(xy 2 + x 2 y - xy) - y 2 z (y - l)(xz 2 + x 2 z - xz) 
-x y z (l + x )(zy 2 + z 2 y - zy) + x l(X 5 )z 2 y 2 

w(x, 2/, z, x , 2/0, zq) = l(X 5 )xy 2 z + (2/0 - 2/o) (xz 2 + x 2 z - xz) - xly (zy 2 + z 2 y - yz) 
-y z 2 (xy 2 + x 2 y - xy) + y (2x y + 2y z + 2x z - 3x - 2y - 3z + 2)xyz 

and l(X 5 ) = Xq + 2/0 + Zq — 1. It is easy to see that if any four nodes do not lie on 
hyperplane then Sj = II3. Hermite interpolation of type total degree is affinely invariant 
in the sense that if the interpolation is singular or regular. Hence if the given five nodes 
are in general position, that is, no four nodes lie on a hyperplane, then any four of them 
can be transformed into Xi, X2, X3, X4. This example implies that uniform Hermite 
interpolation of type total degree on 5 nodes up to order 1 in R 3 is almost regular. 
Example 2. Consider the case of d — 3, m — 6 and p = {3, 3, 3, 3, 3, 3}. Take 

X, = (0, 0, 0), X 2 = (1, 0, 0), X 3 = (0, 1, 0), X 4 = (0, 0, 1), X 5 = (1, 1, 1), X 6 = (2, 1, 1) 

and % = {X 1 ,X 2 ,X 3 ,X 4 ,X 5 ,X 6 }. With the help of Maple, we have S r = U 3 7 . Thus 
uniform Hermite interpolation of type total degree on 6 points up to order 3 in IR 3 is 
almost regular. 
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3 Singular Interpolation Scheme 



In this section, we will investigate Hermite interpolation of type total degree which 
is singular. Our results covers those appeared in [llj , but the method employed here is 
simpler. 

In this section, Eq. (J2J) is always assumed to hold. In this case, the interpolation 
space and the set of functionals to be interpolated are affinely invariant. 

The following theorem and corollary will give an evaluation of n in Eq. (J2]). 

Theorem 4. Given 2£ = {Xi,X 2 , . . . ,X m } and p = {p±,P2, ■ ■ ■ ,p m }, if there exists 
an n such that (j2T,p,n^) correct, then the following inequality holds: 

n + d\ (Vm + d 

where 1 < d < d and the right side denotes the sum of any d + 1 terms. If m < d + 1, 
we assume 

i=i x ' i=i x 



Proof. We only give the proof of m > d+1. The proof of m < d + 1 is similar and 
easier. Note that Eq. (T5]) holds since p, 11^) is correct. Thus inequality (IE]) trivially 
holds for d = d. 

Consider the case of d < d. Suppose X qi ,X q2 , . . . , X q . j are d+1 arbitrary nodes. 

Then there exist d — d linearly independent linear polynomial k(X), i — d + 1, 2, . . . , d 
vanishing on these d + 1 nodes. Assume li(X),i = l,...,d are d linear polynomial 
such that all k(X), i = 1,2, ... ,d are linearly independent. Take the following affine 
transformation 

T: Vi = k{X),i = l,2,...,d (9) 

Let Y{ = T(Xi),i = 1,2, ... ,m and W = T(JT). Thus under the new coordinate 
system, the last d — d coordinates of Y q ., i = 1, 2, . . . , d + 1 are zero. 

Hermite interpolation of type total degree is affinely invariant in the sense that if 
the interpolation is singular or regular. Hence (^, p, n^) is also correct. Thus for any 
given {c q ^a, < |ot| < p q , 1 < q < m} there is a unique /ellf satisfying 

Qcti+ct2+...+a d 

dy^...dyr f{Yq) = Cq ' 0n l - q ~ m > °^H^ ( 10 ) 
Specially, we have 

tt^ ^-ixjf(Yq) = c q , a , l<q<m, <\a\=a 1 + a 2 + ... + a i <p q 

dy^.-.dy/ 
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which means 



( 



n + d 
d 



) 




The proof is completed. 



□ 



For convenience, we always order < p\ < p 2 < • • • < p m in what follows. 



Corollary 1. If set d = 1, then 



n + l>p m + p m „i + 2, or n>p m + p m _i + l. 



(11) 



The following theorem can be use to judge whether the interpolation scheme is 
singular for small m. 

Theorem 5. Given 2£ = {X\,X 2 , . . . ,X m } and p = {pi,p 2 , ■ ■ ■ ,Pm}, if 



then Hermite interpolation of type total degree is singular. Here the numbers p q and n 
are assumed to satisfy Eq. ([2]). 

Proof. We only need to find a polynomial satisfying the homogenous interpolation con- 
dition, which can be done by giving an algorithm for its construction. 

Step 1. Set f{xi,x 2 , ...,x d ) = l and p = {pi,p 2 , ■ ■ ■ ,p m } ■= {pi + l,p 2 + 
l,...,p m + l}. 

Step 2. If the number of the nonzero in p is no more than d, let l(x±, x 2 , ■ ■ ■ , Xd) 

be a linear polynomial vanishing on X\, X 2 , . . . , X m . Take p max = max{p 1 ,p 2 , . . . , p m } 

and set / = / • /P max . Otherwise, go to step 3. 

Step 3. Suppose p^yp^, . . . ,p id are d largest numbers in p. Clearly, there 
must exist at least one linear polynomial vanishing at any d points. De- 
note by 1^2. ..i d the linear polynomial vanishing on X ii: X i2 , . . . , X id . Set 
f = f ■ lhi 2 ...i d - Let pi. = pi 3 = 1,2, ...,d and pi = pi if % e 

{1,2,..., m}/ i 2 ,..., id}- Go to Step 2. 

We want to show that the polynomial / constructed by this algorithm satisfied our 
requirement. Firstly, to this end we need to show that the algorithm does eventually 
terminate. Denote |p| = Y^iLxPi- The key observation is that |p| will be dropped by d 
after step 3. Hence the algorithm will terminate since at the beginning |p| < nd. 

Next, it is easy to know that this kind of polynomial / constructed by the algorithm 
satisfies the homogenous interpolation conditions by Theorem [3j 

Finally, we need to show that deg(/) < n. Assume that Step 3 has been run t times 
totally. Clearly, t is no more than n according to the assumption of the theorem. To 
complete the proof, now we estimate the degree of the polynomial /. Suppose we are 
in a situation to run the final step. That is, there are only no more than d nonzero 
numbers in p. We consider the following three cases. 



Pi + p 2 + . . . + p m + m < nd, 



(12) 
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1) If all the numbers in p are zeros, then the degree of / is t which is not larger than 



n. 



2) The largest number in p equals to 1, that is, ]5 max = 1. In this case, we have 
t < n — 1, which will leads to deg(/) < n. 

3) If p max > 1, we will show that the degree of / is no more than k := max{pi + l,p 2 + 
1, . . . ,p m + 1}. Clearly, it is enough to show that p max + t = k. To this purpose, 
denote by p^ the set p after running the third Step j times. For convenience, 
we also write {pi + l,p 2 + 1, ■ ■ ■ ,p m + 1} as p^ ^. Based on these notation we 
have that p max is the maximum number in pW. Furthermore, there are at most 
d numbers in pW different from zero. Thus we can deduce that p max + 1 must 
be the maximum number in p^ t_1 \ If it is not the case, then there must exist 
d numbers in p( t_1 ) are larger than or equals to p max , which will leads to more 
than d numbers in pW different from zero and furtherly contradict the conclusion 
above. By a similar discussion, we finally derive that p max + t is the maximum 
number in p(°) which implies that p max + t = k. 

By collecting above discussion, we get that / satisfies the interpolation condition and its 
degree is no more than n. Thus the interpolation scheme is singular, which completes 
the proof. □ 

Corollary 2. Given = {Xi,X 2 , . . . ,X m } andp = {pi,p 2 , ■ ■ ■ ,Pm}, if Eq. (fl2"|) holds 
then there exists a polynomial of degree < n, together with all of its partial derivatives 
of order up to Pi, vanishing at X 4 for all i. 

Theorem [5] is sharp for small m. By corollary [1] and theorem [5l we have 

Theorem 6. All Hermite interpolation of type total degree are singular in W d with 
d > 2 if the number m of nodes satisfies 2 < m < d + 1, except for Lagrange case. 



Proof. We only need to show that Eq. (U21) holds in this case. By Eq. (jTTj) in corollary 
[TJ we have n > p m + p m -\ + 1, which will lead to 

PI+P2 + ■ ■ ■ + Pm + m < px + P2 + ■ ■ ■ + Pm + d + 1 

< + p m + d + 1 

< d(p m +p m -l + 1) 

< nd. 

This completes the proof by Theorem [5j □ 



This result appear in H[ and proved for d — 2,m — d + 1 in [Uj and for d > 2 in 



. The following result is more sharp. 



Corollary 3. All Hermite interpolation of type total degree are singular in M d with 
d>k(l + ^) ifm<d + k. 
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Proof. We only need to show that Eq. (fT2"j) holds in this case. Again by Eq. (fTTj) in 
Corollary [p, we have n > p m + p m -i + 1- Thus, 

P1+P2 + ■ ■■ + Pm + m < pi + p 2 + ■ ■ ■ + p m + d+ k 

< (d + k -l)p m -i+Pm + d + k 

< dp m ^i + kp m + d + k 
= dp m _ x + k(p m + l)+d 

= dp m _i + kp m (l H ) + d 

Pm 

< dp m _i + dp m + d 

< nd. 

This completes the proof. □ 
For true Hermite interpolation, p m > 1 which means 1 + — < 2. Thus we have 

Corollary 4. All Hermite interpolation of total degree are singular in M. d with d > 2k 
if the number m of nodes satisfies m < d + k. 

From corollary HI we know that all Hermite interpolation of total degree are singular 
with d > 4 if the number of nodes satisfies m < d + 2. And also from corollary [31 all 
Hermite interpolation of total degree are singular for d = 3 and m = 5 if p m > 2. We 
claim that this result is very sharp because the corresponding Hermite interpolation 
with d — 3,m = 5 and p\ = p 2 = P3 = p$ = Pts = 1 is almost regular. The regularity 



was proved by the method of determinant in [llj and also be shown by example 2 in 
section 2. Besides, if p m < 1 and p\ = 0, Eq. ([2]) never holds. Therefore, for d = 3, 
only one case can produce regular interpolation scheme. 

Now let us consider the case of d = 2 and m = d + 2. It is well known, interpolating 
the value of a function and all of its partial derivatives of order up to p at each of the 
three vertices of a triangle as well as the value of the function and all of its derivatives 
of order up to p + 1/p— 1 at a fourth point lying anywhere in the interior of the triangle 
by polynomials from ^.^+2/^^+1 is regular. We will prove that in all other cases, the 
corresponding Hermite interpolation scheme is singular. This can be done by Lemmas 
SandU 

Lemma 2. All Hermite interpolations of type total degree on 4 points in M 2 are singular 
if0<p 1 <p 2 <p 3 < Pa, P4 > P2 and p 3 > pi . 

Proof. Suppose the interpolation problem is regular, that is, there is a unique f(xi, x 2 ) G 
satisfying 

^ ' " f{X q ) = c q>a , 1 < q < 4, < |a| < p q . (13) 



It follows from corollary [T] that n > p 3 + p 4 + 1 
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Denote by = the line passing through Xi and Xj. Consider the following 
polynomial 

f _ 7Pi+l;P3+l/ max {P4-P3,P2-pi} 
J — '12 '34 '24 

Clearly, / satisfies the homogenous interpolation conditions in Eq. f|T3|) and has the 
degree of 

(pi + 1) + (p 3 + 1) + max{p 4 -p 3 ,p 2 -Pi} = m&x{p 1 +p 4 ,p 2 +p 3 } + 2. 
To complete the proof, it remains to prove n is not less than the degree of /, that is, 

n > max{pi + p 4 , p 2 + p 3 } + 2. 
According to the discussion above, it is enough to show 

Pi + P3 + 1 > max{pi + p 4 , p 2 + p 3 } + 2 
which is equivalent to 

(Pa - P2) + (Ps ~ Pi) > KP4-P2) - (Ps-Pi)l +2. (14) 

Inequality pit! will be satisfied if p 4 > P2 and j»3 > pi. This completes the proof. □ 

Lemma 3. All Hermite interpolation of total degree with d = 2 and m = 4 are singular 
if the orders satisfy one of the two conditions 

0- Pi=P2=Pz andp 4 >p! + 2, 
ii). p 2 = p 3 = Pa and Pi>pi + 2. 

Proof. We only give a proof for case i). In this case we have 

Pi + P2 + P3 + Pa + 4 = 3p 3 + p 4 + 4 

< 2p 3 + 2p 4 + 2 
= 2{p 3 +p 4 + 1) < nd 

which completes the proof by theorem [51 □ 
It is well known that uniform Hermite interpolation of type total degree never hap- 



pen because Eq. fl2]) in this case does not hold, see [111 ]. Thus for m = d + 2, we 
have 

Theorem 7. Consider the problem of Hermite interpolation of type total degree on 
m = d + 2 nodes in M. d . Then 

• For d = 2, if pi = p 2 = P 3 ,Pa = P3 + 1 or pi = p 2 — l,p 2 = p 3 = p 4; it is almost 
regular. 

• For d = 3, if pi = p 2 = p 3 = p 4 = p 5 = 1 and n = 3, it is almost regular. 
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• Otherwise, it is singular. 

Let us consider the case of m = d + 3. From corollary HJ all Hermite interpolation of 
total degree are singular in one of the following cases: i) d > 6, m = d + 3 and p m > 1; 
ii) d > 5, m = d + 3 and p m > 2; hi) d > 4, m = d + 3 and p m > 3. 

For the case of d = 5, m = 8 and p m = 1, it is easy to check Eq. ([2} never holds. 

Let us turn to the case of d = 4, m = 7 and p m < 2. Eq. (j2J) holds only if i) 

= Vi — 2, p« = 0, 1 < i < 5 and n — 3; ii) pe = pt = l,pi — 0, 1 < z < 5 and n = 2. 
For these two interpolation schemes, we need the following result from |l3[ : 

Theorem 8. Multivariate Hermite interpolation of type total degree (8) in R d with at 
most d + 1 nodes having p q > 1 is regular a. e. if and only if 

Pq+p r <n 

for 1 < q,r < m,q ^ r. 

Obviously the interpolation problems considered above are singular by the theorem 
above. 

Consider the case of d = 3 and m = 6. According to Corollary (TJ n > +ps + 1. 
Thus if pq — 1 = j?5 = p4 = ps = P2 = pi = p, then n > 2p + 2 and 




= -_(2p + 2)(p + 2)(p + l) <0 
o 

which implies that Eq. (j2J) never holds. Else if P6 > P5 > Pi, then 

Pi + £>2 + • • • + Pe + 6 < 5p 5 + p 6 + 5 

< 3p 5 + 3p 6 + 3 

< 3(p 5 +P6 + l) <3n. 

Thus it follows from Theorem [5] that the Hermite interpolation problem is also singular 
in this case. 

Otherwise, p 5 = p 6 . In this case 




= --(p 6 -3)(p 6 + 2)(p 6 + l). 

Thus Eq. d2J) does not hold for p 5 = p 6 > 3. Moreover, by a careful check and 
computation, Eq. (|2J) also does not hold for p 5 = p 6 < 3. As a result, Eq. (j2J) only hold 
for pj = 3 and n = 7. By Example 2 in section 2, the uniform Hermite interpolation 
problem is almost regular. 
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Finally, we consider the case of d — 2 and m — 5. In this case, we claim that 
n < 2p§ + 2 if the Hermite interpolation problem is regular. In fact, for any 5 points 
in the plane, there must exist a non-trival quadratic Q(x,y) vanishing at these points. 
Let P(x,y) = [Q(x, y)] P5+1 . Then P, together with all of its partial derivatives of order 
up to p 5 , vanish at these points. 

Lemma 4. Given = {X\, X 2 , . . . , X 5 } C M 2 and p = {pi,p2, ■ ■ ■ ,Ps} with pi < 
< p 3 < p 4 < p 5; if Eq. 02]) holds and 

p 2 +p 3 + 2 < pi + p 5 (15) 

then Hermite interpolation of type total degree is singular. 

Proof. As discussed above, there exists a non-trival quadratic Q(x, y) vanishing at these 
five points. Since 

(?2 - Pi - 1) + (P3 - Pi - 1) + (p 4 - Pi - 1) + (PS ~ Pi - 1) + 4 

= P2 + P3 + Pa + P5 - 4pi < 2(p 4 + p 5 - 2pi - 1), 

due to Corollary [2] there exists a polynomial /(x, of degree < P4+P5— 2p% — 1, together 
with all of its partial derivatives of order up to Pi —p\ — 1 (if negative, no interpolation 
happens), vanishing at Xj for all 2 < % < 5. Let P(x,y) = [Q(x, y)] Pl+1 ■ f(x,y). Thus 
P, together with all of its partial derivatives of order up to pt, vanish at X{ for all i. 
It is easy to get that the degree of P is no more than P4 +P5 + 1. This completes the 
proof. □ 

Lemma 5 ([Hj]). ^ e uniform Hermite interpolation of type total degree on 5 nodes in 
R 2 is singular. 

Lemma 6. If pi < p 2 + 1 = P3 = Pa = Ps and Eq. (J2]) holds, then the Hermite 
interpolation of type total degree is singular. 

Proof. Suppose the interpolation problem is regular. Then, p 4 + p 5 + 1 < n < 2p 5 + 2, 
that is, n — 2p 5 + 1. However, 

Pi + 2\ A f Pi + 2\ fn + 2 



E , >E 



, 2 ^ \ 2 } V 2 

i=l v 7 i=2 v 7 v 

which contradicts Eq. ([2]). □ 

Lemma 7. // pi < £>2 — Pz — Pa — Pa ~ 1 one? .Eg. (J2J) holds, then the Hermite 
interpolation of type total degree is singular. 

Proof. Suppose the interpolation problem is regular. Then we have 

Pa + P5 + 1 < n < 2f>5 + 2 

and 

Pi + 2\ A f P i + 2\ (2p A + 2 + 2 



Err >E 



2 / ^ V 2 / V 2 

8=1 x 7 1=2 
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Thus, n = 2p 5 + 1. Let Q(x,y) be the quadratic polynomial vanishing at these 5 
points. Again by the way of Lemma HI we can get a polynomial of degree no more than 
n — (2pi+2), together with all of its partial derivatives of order up topj— pi — 1, vanishing 
at Xi for 2 < i < 5. Thus [Q(x, y)] Pl+1 ■ f satisfies the homogenous interpolation 
condition. It is easy to check that 

2{ Pl + 1) + n - {2p x + 2)=n 

which completes the proof. □ 
By collecting the discussion above, we have 

Theorem 9. All Hermite interpolation of total degree onm = d + 3 points are singular 
in M. d except for the case of d = 3, n = 7 and Pi = 3 for i = 1,2, ... ,6. 

Combing Corollary H] and Theorems l6]P71l9l we have proved 

Theorem 10. Consider the problem of Hermite interpolation of type total degree on 
m = d + k nodes in lR d . Then 

1. For k < 1, it is singular. 

2. For k = 2, if d = 2, and pi = p 2 = Pz.Pa = P3 + 1 or pi — p 2 - 1,P2 — Pz — Pa, 
it is almost regular; else if d = 3, and p^ — l,i — 1,2, 3,4, 5, n = 3, it is almost 
regular; otherwise it is singular. 

3. For k = 3, if d = 3 and pi = 3,(i = 1,2, ... ,6), n = 7, it is almost regular; 
otherwise it is singular. 

4- If d> 2k, it is singular. 

4 Conclusion 

In this paper, we consider the singular problem of multivariate Hermite interpolation 
of total degree. We make a detailed investigation for Hermite interpolation problem of 
type total degree on m = d + k nodes in M d . Our results imply that the interpolation 
problem in ~R d is singular on m = d + k nodes for d > 2k. For k < 3, it is shown that 
only very few cases can produce regular interpolation. The method developed in this 
paper can deal with the case of small k compared to d. For bigger k, Theorem [5] is not 
very sharp. 
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